Generator Coordinate Calculations for Breathing-Mode Giant 
Monopole Resonance in the Relativistic Mean-Field Theory* 



M.V. STOITSOV a '0, P. RING a and M.M. SHARMA 6 

a Physik Department, Technische Universitat Miinchen, D- 85748 G arching, Germany 
b Max Planck Institut fur Astrophysik, Karl-Schwarzschild-Strasse 1, D-85740 Garching, Germany 



Abstract 



The breathing-mode giant monopole resonance (GMR) is studied within the 
framework of the relativistic mean-field theory using the Generator Coordi- 
nate Method (GCM). The constrained incompressibility and the excitation 
energy of isoscalar giant monopole states are obtained for finite nuclei with 
various sets of Lagrangian parameters. A comparison is made with the re- 
sults of nonrelativistic constrained Skyrme Hartree-Fock calculations and with 
those from Skyrme RPA calculations. In the RMF theory the GCM calcula- 
tions give a transition density for the breathing mode, which resembles much 
that obtained from the Skyrme HF+RPA approach and also that from the 
scaling mode of the GMR. From the systematic study of the breathing-mode 
as a function of the incompressibility in GCM, it is shown that the GCM 
succeeds in describing the GMR energies in nuclei and that the empirical 
breathing-mode energies of heavy nuclei can be reproduced by forces with an 
incompressibility close to K = 300 MeV in the RMF theory. 
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I. INTRODUCTION 



The nuclear matter incompressibility signifies an important and cardinal point on the 
equation of state (EOS). The behaviour of the nuclear matter at the saturation point is 
relevant not only to the property of finite nuclei, but also to astrophysical phenomena such 
as supernovae explosion and neutron stars. The breathing-mode giant monopole resonance 
(GMR), whereby nuclei undergo radial density oscillations, provides a source of extracting 
the dynamical behaviour, i.e., the compression properties of nuclei and nuclear matter 0. In 
addition to the GMR excitation mode which represents a small-amplitude collective motion, 
the intermediate energy heavy-ion collisions |§, on the other hand, strive to map out the 
EOS of the nuclear matter for densities higher than the saturation density. This is also 
expected to constrain the incompressibility at the saturation point. However, owing to the 
complex interplay of many degrees of freedom in the heavy-ion collision, it has not yet been 
possible to gain much insight into the behaviour of the EOS. For properties around the 
saturation point, the GMR remains an important object of investigations. 

The GMR has been measured almost all over the periodic table Some time ago, the 
GMR energy was obtained || with a considerable precision in a set of medium heavy Sn 
and Sm nuclei. Attempts were made to extract the nuclear matter incompressibility from 
such precision measurements. An earlier analysis based upon a leptodermous expansion of 
finite nuclear incompressibility into various finite-size components, led to the nuclear matter 
incompressibility of 300 ±25 MeV |J. This analysis, which took into account the correlation 
of the Coulomb term involving the third derivative of the EOS, was based upon the sys- 
tematics from the density-dependent Skyrme interactions. In a real case, however, Skyrme 
forces might not be reliable for this purpose. Circumventing this constraint based upon the 
Skyrme interactions, it was found that error bars on the nuclear matter incompressibility in- 
creased by more than 50% and the value itself was obtained at slightly higher than 300 MeV 
0. More recently an analysis of experimental data including deformed nuclei and of data 
from many laboratories was attempted [§]. However, this analysis, which comprises data of 



2 



various origins, was not conclusive on the extraction of the nuclear matter incompressibility. 
A detailed and critical analysis of empirical breathing-mode GMR data is in progress. 

Theoretically, the incompressibility has been obtained using the density-dependent inter- 
actions ||. The deductions base themselves upon an interpolation between various Skyrme 
and Gogny forces for the GMR energies obtained from self-consistent HF+RPA calculations. 
These calculations were a major effort with a view to explaining the breathing-mode energies 
in finite nuclei in a microscopic approach. This approach, however, succeeded in reproduc- 
ing the GMR energy of only 208 Pb within the interpolation scheme. The GMR energies of 
90 Zr were overestimated by 1-2 MeV. This fact has been corroborated by the calculations 
within the RPA sum-rule approach using various Skyrme interactions ||. The calculations 
indeed reproduce the GMR energy of 208 Pb using Skyrme force SkM*. The GMR energies of 
medium-heavy nuclei such as 90 Zr, Sn and Sm isotopes could not, however, be reproduced 
within the Skyrme forces. The Skyrme interaction SkM* has been used extensively to calcu- 
late the properties of giant resonances jnj. It reproduces the empirical excitation energies 
of giant quadrupole resonance (GQR) very well. The appropriate effective mass of this force 
helps to achieve the required GQR energies. The force SkM*, however, reproduces the GMR 
energies of only 208 Pb well. This is due to a simple relationship of the surface incompress- 
ibility to the bulk incompressibility for the Skyrme type of forces, that for a given force the 
surface incompressibility has about the same value as the bulk incompressibility (Tl|. This 



relationship has essentially been at the root of the problems in describing adequately the 
mass dependence of the GMR energies in the Skyrme ansatz. 

Relativistic mean-field (RMF) theory Jl2| has in the last years been found to be espe- 
cially appealing in describing the ground-state properties of nuclei at and far away from the 
stability line [0,^]. The long-standing problem of the kink in isotope shifts in Pb nuclei, 
which could not be described with the Skyrme forces including all possible correlations, has 
been successfully solved in the RMF theory ||15|| . The theory has subsequently also been 
able to provide a good description of the binding energies and deformations of nuclei close 
to neutron drip-line ||16|| . Shell effects arising from the Dirac structure of the spin-orbit in- 



teraction in the RMF theory manifest in the behaviour of the binding energies. The strong 
shell-effects arising from the RMF theory are corroborated by the finite-range droplet model 
(FRDM) |T7|] and are in contrast with those from the Skyrme theory [|]1||T9|]. Thus, the RMF 



theory has achieved a considerable success in describing many aspects of the ground-state 
properties of nuclei. 

The dynamical aspects within the RMF theory have remained largely unexplored. A first 
attempt was made to obtain the breathing-mode energies and incompressibilities within the 
RMF theory using the linear Walecka model in constrained calculations [pOR . Such calcu- 
lations were further extended to light nuclei and anharmonicities in the breathing-mode 



oscillations were indicated |^TJ. The relationship of the GMR energies to the incompress- 
ibility of nuclear matter is, however, not yet known for the RMF theory. On the contrary, 
in the Skyrme approach, the relationship between the GMR energies and the incompress- 
ibility has been studied extensively (see e.g., refs. PJlUflfljlOll) and has been found to be 
straightforward. An exercise to understand this relationship in the RMF theory has recently 
been undertaken [^] employing relativistic constrained calculations within the mean field. 
Another approach which has received considerable attention as a useful tool to study prop- 
erties of excited stated in nuclei is the generator coordinate method (GCM) J23J. It has been 
applied amongst others also for the breathing mode |24Tf26f. This has been attempted in 



the non-relativistic theories with a view to taking into account the relevant correlations in 
the nuclei. In this paper, we investigate the GCM for the first time in the RMF theory and 
focus upon the structure and properties of the breathing-mode GMR using the method of 
generator coordinates. A comparison of the properties of the GMR will be made with those 
from the Skyrme ansatz. 

The paper is organized in the following way: In Section II we provide the theoretical 
framework of the RMF theory. The details on the Generator Coordinate Method in the RMF 
theory are presented in Section III. The problem of the breathing mode GMR is discussed in 
Section IV. In Section V we discuss the results obtained in this framework. The last section 
contains a summary and conclusions. 



II. RELATIVISTIC MEAN-FIELD THEORY 



We start from Relativistic Mean Field theory [12] which treats the nucleons as Dirac 
spinors ip interacting by the exchange of several mesons: scalar cr-meson that produces 
a strong attraction, isoscalar vector w-meson that causes a strong repulsion, isovector p- 
meson required to generate the required isospin asymmetry and photon that produces the 
electromagnetic interaction. The model Lagrangian density is: 

C = - M}tp + \d^ad p a - U(a) - g^crip 

-in^fya, + \mluj^ - g^i^u^ 

where U(a) is the non-linear scalar self-interaction with the cubic and quartic terms required 
for appropriate surface properties 



U{a) = \rnla 2 + l -g 2 <y* + \g z a\ (2) 

M, 77V, Tn w and m p are the nucleon, the ex-, the lu-, and the p-meson masses, respectively, 
and g a , g u , g p and e 2 /47r=l/137 are the coupling constants for the a-, the o>-, the p-mesons 
and for the photon. The field tensors for the vector mesons are: 

QT = d»uo v - d u ou>*, (3) 

BT = - d v F - g p {^ x p) (4) 
and for the electromagnetic field 

The associated Hamiltonian operator H is then obtained using the well known canonical 
quantization procedure based on the anti-commutator (for the fermions) and the commutator 
(for the mesons) relations |28|j29| . 



Within the relativistic mean-field (RMF) approximation the A independent nucleons with 
single-particle spinors ipi, (i — 1, 2, A), are assumed to form a single Slater determinant 
$ and to move independently in the meson fields. In the particular case of spherical nuclei, 
symmetries simplify the calculations considerably and only the time-like components u°(r), 
p°(r) and A°(r) of the u-, the p- and the electromagnetic fields survive. When describing 
groundstate properties of nuclei, one looks for static field solutions 0(r) = a(r), u°(r), p°(r) 
and A°(r) that satisfy the Klein-Gordon equation 



( 92 29 2\ j., , 



(6) 



where are the meson masses for <f> = a, u, p and is zero for the photon. The source 
terms 



—g a p s (r) — g 2 (J 2 {r) — g 3 cr 3 (r) for the a-field 

Q^Pvij) for the tu-field 

9pPz( r ) f° r the p- field 

epp(r) for the Coulomb field, 



(7) 



depend on the spherical densities 

A 



i=l 
A 

Pv(r) =J2^l( r )M r 

i=l 
A 



(8) 



i=l 

pp(t) = £v4(r)^^W) 



i=i 

where, in the no-sea approximation, the summation runs over all occupied states in the 
Slater determinant $. The solution of Eq. can, in principle, be expressed in terms of 
Green's functions, i.e., 



00 

G^(r,r')^(r')r' 2 cir', (9) 







where, for the massive fields, 



G ( r r ') = 11 ( e -m,\r-r'\ _ e -m,\r+r'\\ (1Q) 

Im^ rr v 1 

and for the Coulomb field, 

{1/r for r > r' 
(11) 
1/r' for r < r'. 

The total ground-state energy of spherical nuclei can be expressed, in the center-of-mass 
frame, as a functional of the baryon spinors {ipi} 

Ermf[H = <®\H\$> (12) 
where the Hamiltonian density 

n RMF {r) = r(r) +Mp s (r) + \g aPs {r)a{r) + H|^a 3 (r) + \g^\r)} 

(13) 

+ \g uPv {r)uj\r) + \g P Pz{r)p\r) + \ep p {r)A\r) 
depends only on the baryon field since the 'kinetic' energy density 

r(r) ee $>J(r){-iaV}^(r), (14) 

i=i 

and the spherical densities and therefore the mesonic fields (|9]) are all expressed in terms 
of the Dirac spinors {ipi}- In the RMF approach Fock terms in Eq.flTB]) are neglected. 

Taking the variation of Eq.(^) with respect to ip\ one obtains the stationary Dirac 
equation with the single-particle energies as eigenvalues, 

h D ijji(r) = e f V>i(r), (15) 

where 

h D = -iaV + (3(M + g a a(r))+ giJ LU (r)+g p T 3 p (r) + e£^A (r). (16) 

Solving this equation self-consistently (the mesonic fields depend on the baryon solution 
according to Eq.(|9])) one obtains the nuclear ground state 3>o m terms of the solutions {ipi}. 



7 



III. RELATIVISTIC GENERATOR COORDINATE METHOD 



The GCM has been used extensively within the non-relativistic approaches to obtain 



we present a relativistic extention of the generator coordinate method (GCM), which is 
based upon a trial A-particle wave function ansatz ^gcm written in the form of a linear 
combination: 



where the generating function <&(q) = $(ri, . . . ta] q) is chosen to be a Slater determinant 
built upon single-particle spinors ipi(r,q), {i = 1,2, A), depending on the generator 
coordinate q. It is obvious that in this case the wave function of the system (|17D, being 
a superposition of Slater determinants goes beyond the limits of the RMF approach. 

The so-called 'weight', or 'generator' function F{q) is determined after varying with respect 
to ^(q) the energy of the system 



the ground-state and excited states of nuclei P0fl . Using the Skyrme forces, the GCM was 
applied to study the giant resonances E|. Recently, the GCM has also been employed 
to investigate the effect of correlations on the ground-state properties of nuclei pi] . Here 




(17) 



E[F\ 



< ^gcm\H\^gcm > 

< ^ GCM | ^ GCM > 



(18) 



This leads to the Hill- Wheeler integral equation for the weight function: 




(19) 



where 



H(q,q') =< <%)|tf|<%')> 



(20) 



and 



X(q,q') =< <%)!<%')> 



(21) 



are the energy and the norm overlap kernels, respectively. 



S 



A straightforward calculation shows that with the Hamiltonian H associated with our 
model Lagrangian (P one obtains 



H(q,q') = <$(q)\HMq')> 



(22) 



where Af(q, q') is the overlap kernel (^) and TC(r; q, q') is the overlap energy-density kernel: 

H(r;q,q') = r(r; q, q') + Mp s (r; q, q') + ^g a p s (r; q, q')a(r; q, q') 

+Hl92(T 3 (r; q, q') + \gz(y A {r; q, q')} + \g^p v {r\ q, q')uj°(r; q, q') (23) 
+l9 P P3(r; q, q')p°(r; q, q') + \ep p (r; q, q')A°(r; q, q'). 

In this equation the 'kinetic' energy density is defined by the spinors {ipi(r; q)} as 



r(r;q,q')= £ ^ ^J(r; g){-iaV}^(r; ?')■ 

*j=i 



Similarly the other densities entering Eq.(|23"D are: 



(24) 



p s (r;q,q') = J2 N ji* Q^jfa l') 

A 

pv(r-,q,q') = J2 N 7i~ Wto ?)^-( r ; i') 



P3(r;q,q') = J2 fyi* V>l( r ; 4 



i,j=l 

A 



p p (r;q,q>) = £ TV/ ^(r; g)^^-^-(r; g'). 



They appear as source terms 



s^(r;q } q' 



(25) 



-9aPs(r; q, q') - g-io 1 - S^o" 3 for the cr-field 

9uPv( r ] q, q') for the w-field 

g p Pz{r \ q, q') for the p-field 

ep p (r; q, q') for the Coulomb field, 

in Klein-Gordon equations of the type (|6|) whose solution determines the fields </>(r; q, q') 
cr(r; q, q'), u)°(r; q, q'), p°(r; q, q') and A°(r; q, q') entering Eq . (^3|) as 



(26) 



[r;q,q)= I G^(r, r')^(r'; q, q'y 2 dr', 
/o 



(27) 



the Green functions are defined as before by Eqs. (0) and (pTTj) . 

In the above equations the sums run over all occupied single-particle states and N^ 1 are 
the elements of the matrix N~ l (q, q') where 

N tJ (q,q') = J dh4{v-q)^{v-q'). (28) 

The determinant of N(q, q') simply gives the overlap kernel (|2~T|) 

Af(q,q')=det{N(q,q')}. (29) 



Thus, having determined the integral kernels (|22|) and (p9|), the associated Hill- Wheeler 
integral equation ([19]) has to be solved in order to determine the nuclear ground- and n- 
excited states through its eigen solutions {E , ^(q)} and {E n , J r n {q)} 1 respectively 

IV. ISOSCALAR GIANT MONOPOLE RESONANCE 

The constrained Hartree-Fock calculations have been a usual method to obtain descrip- 
tion of the excited states in nuclei. An extension of this method in the framework of the 
RMF theory has been made recently, where the breathing-mode GMR in finite nuclei has 



been obtained in the constrained calculations [25]. We extend some of the discussion here 
for the sake of clarity. In order to analyze the isoscalar GMR we perform also constrained 
RMF calculations, where the Dirac equation (see Eq. (|15|)) 

(h D - gr 2 ) ipi(x) = e#l)i(x), (30) 

is solved at different values of the Lagrange multiplier q which are associated with values of 
the nuclear rms radius 

R = {jjr 2 p v (r;q)d 3 r} 1/2 , (31) 

where 

A 



1=1 



Pv(r; q) = J2 ^ ( r > <?)?Mr; q) (32) 
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is the baryon local density determined by the solution {ipi(r; q)}. According to Eq.flTJ), the 
total energy of the constrained system 

E RM F{q) = ERM F [ipi(q)], (33) 

is a function of q (or the nuclear rms radius R). It has a minimum, the ground-state 
energy E RMF = Ermf(0), at g = corresponding to the ground-state rms radius i?o- The 
curvature of this function around the equilibrium point Rq defines the so-called constrained 
incompressibility coefficient of the finite nucleus 

Kc{A) = A" 1 [r ^^I^ ) ■ ( 34 ) 

The constrained energy ([33]) as a function of q represents the energy surface for the isoscalar 
monopole motion of the nucleus where R changes around its ground-state value Rq. In order 
to derive vibrational excitation energies one needs in addition the inertial parameter for this 
motion. In the non-relativisitic RPA sum-rule approach (SRA) |32| the inertia parameter 
for the GMR is derived as MRq. In this case one obtains the GMR excitation energy E\ as 



Eh 



(35) 



\ MRq ■ 

In order to obtain a description of the GMR in the RMF theory, we consider the Lagrange 
multiplier q entering Eq. (|30D as generator coordinate for the GCM calculations as described 
in Sec.|THl The solution {ipi(r; q)} of Eq.(|30D at different values q then defines the generator 
Slater determinants $(<?) and therefore the integral kernels ( p0|) and (pip . In fact, the 
diagonal part of the energy kernel TC(q,q') coincides with the constrained energy fl3"Bp, i.e., 
E RMF (q) = H(q,q). The off-diagonal elements TC(q,q') contain the information about the 
inertia. We then solve the resulting Hill- Wheeler equation ( |l9l) numerically using the method 



of ref. 25 . 



V. RESULTS 
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A. Relativistic GCM Calculations 



We have performed GCM calculations for four closed-shell nuclei 16 0, 40 Ca, 90 Zr and 
208 Pb with the sets of Lagrangian parameters given in Table|. The Lagrangian parameters 
sets are NL1 J£J, NL-SH [|TJ, NL2 |4|, HS @ and LI @ in the increasing order of 
the nuclear matter incompressibility with K^m =211.7, 355.0, 399.2, 545 and 626.3 MeV, 
respectively. This allows us to examine the dependence of GMR energies on the nuclear 
matter incompressibility AOvm- These sets of parameters have also been employed in our 
earlier constrained RMF calculations [|22[| . The last two sets, HS and LI, correspond to the 
linear model without the self-coupling of the cx-field. In addition, the set LI excludes the 
contribution from the p-field. Among the sets NL1, NL-SH and NL2, which correspond to 
the non-linear model, only the set NL2 has a positive coupling constant in Eq (2). The 
set NL2 has an effective mass m* = 0.67 at the saturation, which is higher than that of NL1 
and NL-SH. Whereas the set NL1 reproduces the ground-state properties of nuclei only close 
to the stability line due to the very large asymmetry energy, the set NL-SH describes also 
nuclei far away from the stability line |14| . The shell effects at the drip-line and deformation 
properties obtained with NL-SH have been found to be in good agreement with the recent 



finite-range droplet model ||17|| . In order to cover the region of the incompressibility about 
300 MeV, we have also constructed a schematic force NL-S1 with K^m — 296 MeV. This 
force describes the ground-state properties of closed-shell nuclei rather well, but has a large 
asymmetry energy of 52 MeV. The use of various parameter sets in our GCM calculations 
allows to study the influence of the model Lagrangian and their properties on the properties 
of the GMR energies. The aim of the present study is also to examine how the GCM works 
in the RMF theory. 

The self-consistent solution of the constrained mean field problem ([30]) diverges for large 
positive and negative values of the constraining parameter q. For large positive values the 
quadratic r-dependence of the constraining operator leads to unbound solutions. For large 
negative values of q, when the nuclear rms radius R decreases, the self-consistent mean-field 
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potential pushes up the single-particle energies and the RMF solution disappears at some 
minimal rms radius R m in- The instability for negative g-values causes no further problems. 
For positive q, however, we find only a very limited range of possible g-values close to the 
ground state. We therefore introduce a cut-off function (1 + exp(r — R cu t/o)~ l ) multiplying 
r 2 in the constrained RMF calculations. The cut-off radius was chosen to be R cut = 2tqA 1 ^ 3 
with ro = 1.2 fm. The diffuseness parameter was set at a = 0.5 fm. In Fig. [I] we show 
the first three excited states for 90 Zr obtained in the relativistic GCM calculations with 
the force NL-SH. The effect of the cut-off function is demonstrated in this figure by two 
constrained RMF calculations, one with (dot-dashed line) and the other without (solid line) 
the cut-off function. The GCM kernel without cutoff goes only up to 4.35 fm on the right 
side of the ground-state as shown by the solid curve. It can be seen that the inclusion of the 
cut-off function enlarges the required space for the GCM calculations without changing the 
behaviour of the integral kernels beyond the one calculated without the cut-off function. 



B. The Correlated Ground State 

Solving the Hill- Wheeler equation (|i~9D we obtain the weight function g = M 1 I 2 Tq 
associated with the ground-state solution To(q). Functions go(q) show the usual bell shape 
with a maximum around the RMF ground-state value q=0. With an increase in mass number 
A the width of go(q) decreases, while its amplitude increases keeping fixed the normalization 
TMT = 1. 

Typical GCM results for energies and rms radii calculated with the set NL1 are given 
in Table |I[ It is worth noting that the GCM ground-state energy is slightly lower than 
the RMF one. This small difference, which is too small to be seen in Fig. 1 contains in 
fact two contributions: (i) the positive zero-point energy of roughly l/2huj in the harmonic 
approximation and (ii) the correlation energy induced by the GCM-correlations lowering the 
mean field energy of the ground state by roughly the same amount. This is an important 
point and reflects the fact that GCM is beyond the RMF approximation. 



13 



There is also no perceptible effect on the rms radii of the nuclear ground state in the 
GCM. The largest difference between the RMF and the GCM ground-state rms radii is seen 
for the nucleus 16 0. It is about 0.0025 fm. Fig.@ shows the RMF and the GCM vector (p v ) 
and scalar (p s ) densities. The RMF and GCM local densities do not differ significantly. For 
heavy nuclei the GCM ground-state local densities are even closer to the uncorrelated RMF 
ones. We can thus conclude that the correlations in the GCM ground-state are small and 
the main purpose of the GCM consideration here is in its possibility to generate nuclear 
excited states. 



C. GMR Excited States 

The first three GMR excited states obtained in the GCM with the parameter set NL1 
are shown in Table [TT]. These states show a clear equidistant spectrum for heavy nuclei. A 
similar behaviour is also apparent from Fig. 1 for 90 Zr too, which has been shown for the 
set NL-SH. For the lighter nuclei, however, there are significant deviations from this type 
of spectrum as can be seen from the energies of the excited states in 40 Ca and 16 0. With 
an increase in mass number, the excitation energies decrease. Here we take the excitation 
energy AE% = Ei—E , which is equivalent to the excitation energy of a collective state in the 
non-relativistic constrained Hartree-Fock approach. The mass dependence AE± = cA^ 1 ^ 3 
of the excitation energy for 208 Pb, the nucleus on which there exists well-measured GMR 
energy, is obtained as c = 69.1, 79.6, 93.6, 104.9, 97.0, and 126.2 MeV for the sets NL1, 
NL-S1, NL-SH, NL2, HS, and LI, respectively. 

In Fig.0 we show the local vector and scalar densities Poo( r ) of the GCM ground state 
and pn(V) of the first excited GMR state for the set NL1. The densities are more extended 
in space in comparison with the ground state ones. Consequently, the rms radii in the first 
exited state are larger than that of the associated ground-state values by about 0.15 fm in 
16 and by only 0.015 fm in 208 Pb. 
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D. Transition Density 



The transition density of the GMR provides the strongest evidence for the radial density 
oscillations in nuclei and hence of the 'breathing' or the compression character of the GMR 
mode. We show in Fig. [3] the vector and scalar transition densities poi( r ) for protons and 
neutrons in 208 Pb obtained in the relativistic GCM calculations for the force NL-SH. The 
transition densities show a change in the density in the bulk at the expense of that in the 
surface. A node at 6 fm is clearly to be seen for protons and at about 6.4 fm for neutrons. 
The existence of a well-defined node in the transition density is the typical behaviour for 
the breathing-mode motion and testifies for the compressional property of the GMR. The 
transition density from RPA calculations for the GMR in 208 Pb with Skyrme force SIII 
was obtained to be very similar to the transition density in Fig. 3. The transition density 
for SIII also showed a node at about 6.2 fm. Both these transition densities, one in the 
RMF theory for NL-SH and the other in the Skyrme approach for SIII resemble much that 
obtained from a simple radial scaling of ground-state density. The point of difference to be 
noted is that in our RMF case, we have obtained the transition density for the GMR in the 
GCM, with some form of a constrained motion. Here we do not observe any conspicuous 
differences between the transition densities of the relativistic GCM and the scaling mode 
in the Skyrme approach. The oscillations in the interior of the nucleus are obviously due 
to the shell effects. The vector transition density shown in the figure conserves the particle 
number. The same can not, however, be said for the scalar transition density, which is albeit 
similar to the vector transition density, but manifests mainly the relativistic effect similar 
to that exhibited by total scalar density. 

We observe that the difference between the scalar and the vector transition densities, 
which is connected with the small components of the Dirac wave functions, arises mainly in 
the interior of the nucleus. In the surface region both densities coincide more or less. 
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E. Constrained Incompressibility of Finite Nuclei 



We now consider the constrained incompressibility Kq{A) as calculated from Eg. ([34]) . 
The results as a function of the nuclear matter incompressibility K^m are shown in Fig. 
^ for a few nuclei. It may be worth mentioning that empirically the GMR has been well 
established only in heavy nuclei, such as 208 Pb and 90 Zr. We have also included the light 
nuclei such as 40 Ca, 16 and 4 He. In the light nuclei it is very uncertain and a full energy- 
weighted sum-rule strength has rarely been observed. Thus, in our case the light nuclei serve 
mostly the purpose of illustration and for the possible anharmonic effects. 

With exceptions for light nuclei, the incompressibility Kc{A) shows a strong dependence 
on the nuclear matter incompressibility Knm- For the linear force HS, Kc{A) shows a slight 
dip from the increasing trend for 208 Pb and 90 Zr, whereas for light nuclei 40 Ca, 16 and 4 He, 
the HS values are even smaller than the NL2 values. The dependence of Kc(A) in the 
Skyrme approach is different, where it increases monotonically with K^m- For, the finite 
nuclear incompressibility receives a sizeable contribution from the surface incompressibility, 
this difference could be explained from the difference in the behaviour of the surface incom- 
pressibility in the two methods. In the Skyrme approach, the surface incompressibility has 
been shown to be Kg ~ —K^m for all standard Skyrme forces. This does not seem to be 
the case for the RMF theory, however, as shown by the HS values. Thus, the surface incom- 
pressibility is not necessarily a straight function of the nuclear matter incompressibility in 



the RMF theory. This point has also been dealt with in ref. |22 



F. Comparison with Nonrelativistic Calculations and Experimental Data 

The excitation energy AEi corresponds to the energy usually obtained from the 
non-relativistic constrained Skyrme Hartree-Fock (SHF) calculations within the sum rule 
approach . In Table [Hj energies AE± are compared with such nonrelativistic constrained 



SHF results obtained with the Skyrme-type forces SkM and SIII. These Skyrme forces have 
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nearly the same nuclear matter incompressibility Knm as do the sets NL1 and NL-SH, 
respectively. It can be seen that the nonrelativistic SHF results differ slightly from the 
values of AE\. This difference in the relativistic GCM excitation energy AE\ from the SHF 
energy is small for heavy nuclei. It, however, increases for lighter nuclei, where the GCM 
shows lower values. 

Fig. [| shows the GCM breathing-mode energy AE\ for various nuclei and parameter sets. 
The energy AEi first increases with Knm from NL1 to NL2 almost linearly for all nuclei. 
For the force HS, which has K^m even larger than that of NL2, the energy shows a decrease 
for all the nuclei, however. This is due to a rather large surface incompressibility which is 
in disproportion to its bulk incompressibility for HS. This reduces the incompressibility of 
the nuclei, as also shown in Fig. 4. For the force LI, AEi shows an increase compared to 
HS. Thus, AE 1 is not related in a simple way to K NM . This reflects the role played by the 
surface component of the compression in the RMF theory. Even for heavy nuclei AE\ does 
not show an overall increasing tendency with K^m- For lighter nuclei this effect is even 
more apparent. 

It is interesting to note that the dip in energy AEi for HS seems to signal the transition 
from nonlinear (NL1, NL-SH, NL2) to linear (HS, LI) models in the Lagrangian (|I]). Even 
with significantly higher nuclear matter incompressibility (Knm= 545 MeV for HS) the 
linear model gives comparable GMR excitation energies (and even lower) in comparison 
with the nonlinear ones (notice that K^m = 399.2 MeV for NL2). 

It is instructive to see that the approximate expression, (|35| ) which is exactly the same 
as in the nonrelativistic sum-rule approach but calculated with the incompressibility Kc(A) 
emerging from the relativistic RMF calculations, gives acceptable results for GMR excitation 
energies. The results from Eq.(p5|) are compared with the GCM and SHF results also in 
Table |IJj. 

In the non-relativistic approach using density-dependent interactions, extensive work was 
carried out to obtain the incompressibility and breathing-mode energy || with HF + RPA 
calculations. The RPA calculations were performed on a set of Skyrme and Gogny interac- 
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tions including the finite-range Gogny force Dl, with an increasing order of incompressibility 
of nuclear matter. This work attempted to reproduce the empirical breathing-mode energies 
on 208 Pb and 90 Zr, where experiments showed the existence of the GMR unambiguously. 
The GMR energy in 208 Pb is rather well-established and lies at 13.7 ± 0.3 MeV. The GMR 
energy in 90 Zr has been measured to be in the range 16.5 - 17.3 MeV by different experi- 
ments. The average value of the energy from different experiment comes at about 17.0 MeV. 
Table [IV] shows the GMR values for 90 Zr and 208 Pb obtained in the RPA calculations |J for 



the forces Dl, Ska and SIII. A comparison of the RPA values with the empirical values in 
Table [IV| shows that the Gogny force Dl reproduces the GMR energy for 208 Pb quite well. 



The force Dl, however, overestimates the GMR energy for 90 Zr by about 1.5 MeV. The 
force Ska, which has incompressibility of nuclear matter at 263 MeV, gives the GMR energy 
for 208 Pb only slightly higher than Dl. The GMR energy with Ska for 90 Zr is, however, 
about 2 MeV higher than the empirical value. Thus, within the non-relativistic approach, 
with Dl one comes very close to reproducing the GMR energy of 208 Pb in the RPA calcu- 
lations. The GMR energy of 90 Zr could not, however, be reproduced by any Skyrme force. 
This has been the scenario within the Skyrme approach, where the conclusions of ref. || 
on the incompressibility hinged very strongly on 208 Pb only. Consequently, a value of the 
incompressibility of nuclear matter of about 210 MeV seem to have been favoured. 

We now compare the empirical values and the RPA results with those in the relativistic 
GCM calculations with the force NL-S1. It may be noted that this force which describes 
the ground-state properties of nuclei only very approximately, and was constructed with a 
view to fill in the gap at about K NM ~ 300 MeV in the the dependence of the incompress- 
ibility on the breathing-mode energy. It has presently only a schematic character. With its 
incompressibility of 296 MeV, the GMR energy for 208 Pb in the GCM has been obtained at 
13.4 MeV. It is very close to the empirical values obtained in many experiments. The GCM 
energy for 90 Zr has been obtained at 17.6 MeV, which is slightly higher than the average 
value of 17.0 MeV but is closer to an earlier empirical result. On the whole, it is within the 
uncertainties of the empirical values. In comparison, the GMR energy from the Gogny force 



Dl in the RPA lies at 18.5 MeV. Systematics of the values for 208 Pb and 90 Zr in Fig. 5 show 
that both the empirical values as shown by the quadrangles are encompassed by the GCM 
calculations curve from K^m = 290 — 310 MeV. The width of the quadrangles signify the 
corresponding experimental uncertainties in the determination of the GMR centroid ener- 
gies. The empirical values themselves have been reproduced by K NM ~ 300 MeV as can be 
seen by intersecting the empirical values at the curves for 208 Pb and 90 Zr. 

VI. CONCLUSIONS AND DISCUSSION 

We have performed a systematic study of the breathing-mode energy and the incom- 
pressibility of finite nuclei with the generator coordinate method in the RMF theory. It has 
been observed that the transition density of the giant monopole mode shows a character 
very similar to that obtained in the Hartree-Fock-RPA approach with density-dependent 
Skyrme interactions. This behaviour is also similar to what one expects in the simple radial 
scaling of the ground-state density. 

Using a set of relativistic mean-field Lagrangian parameters, it has been shown that the 
GCM energies for the realistic forces show an increasing tendency with the nuclear matter 
incompressibility. Only for unrealistic forces such as HS does one observe a decrease in the 
breathing-mode energy and also in the incompressibility of nuclei even when this force has 
a larger K NM . This is due to a very large surface incompressibility of HS. 

The GCM values obtained with the force NL1 are quite lower than the empirical values 
and those with NL-SH are a little higher than the latter. The empirical GMR energies, on 
the other hand, can be well encompassed by the GCM curve from K = 280 — 310 MeV. 
This is corroborated by the GCM values obtained with a rather schematic force having 
an incompressibility K = 296 MeV, where the GCM values for 208 Pb and 90 Zr are very 
close to the corresponding empirical values. Thus, the empirical GMR energies for both 
these nuclei have been clearly bracketed by the GCM calculations in the RMF theory. We 
know of no other theoretical result where the GMR energies for both these nuclei have been 
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reproduced. Our results also bring about severe constraints on the value of the nuclear 
matter incompressibility, the observable which has theoretically been held rather uncertain. 
The GCM results, thus, favour an incompressibility at about 300 MeV. This is in contrast 
with the usual assumption of the incompressibility of about 210 MeV concluded from the 
non-relativistic Skyrme ansatz, where the empirical value for 208 Pb only could be reproduced. 
Our conclusion, on the other hand, is in good agreement with the analysis of the empirical 
breathing-mode energies where the incompressibility of nuclear matter was obtained as 300 
MeV or higher f||7j. This analysis is, however, not yet complete and further work on it is 
in progress. 

Differences in the shell-effects of the RMF theory and the Skyrme approach and their 
implications on the ground-state properties of nuclei such as isotope shifts 1TH] and on 



nuclei at drip lines [|16[ have been discussed earlier. The present work on the breathing- 
mode energies in the GCM has brought about important differences also in the dynamical 
properties of the RMF theory and the Skyrme ansatz. The nearly good reproduction of 
the empirical GMR energies in the relativistic GCM approach has become possible due 
to the ratio of the surface incompressibility to the bulk incompressibility, which has been 



obtained as different from 1 in the RMF forces. This was also demonstrated in ref. |3T) 



using 



various schematic parameter sets that in the RMF theory it is possible to obtain the ratio 
of the surface incompressibility to the bulk incompressibility of up to about 2 or more. For 



the realistic parameter sets NL1 and NL-SH, this ratio has been shown [^2] to be higher 
than 1 (1.58 and 1.72 respectively) in the simple radial scaling of the ground-state density 
in the semi-infinite nuclear matter with the Thomas-Fermi approximation. In the Skyrme 
approach, the ratio of 1 has essentially been at the origin of problems in describing the mass 
dependence of the GMR energies. Of course, there still remain some improvements to be 
made in the ansatz of the RMF theory with a view to describe accurately the ground-state 
energies of nuclei at and far away from the stability line as is the case with the force NL-SH 
as well as the dynamical properties such as the breathing-mode energies in nuclei. 
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TABLES 



TABLE I. Parameter Sets for the Lagrangian (|T|) 



nli m 



NL-SH |4| 



NL2 §4\ 



HS |5| 



li mi 



M a 

m a 

m w 

m p 

9* 

9u 

9p 

92 
.93 



938.0 
492.25 
795.355 
763.0 
10.138 
13.285 
4.975 
-12.172 
-36.265 



939.0 

526.0592 

783.0 

763.0 
10.44355 
12.9451 
4.3828 

- 6.9099 

-15.8337 



938.0 
504.89 
780.0 
763.0 
9.111 
11.493 
5.507 
-2.304 
13.783 



939.0 
520.0 
783.0 
770.0 
10.47 
13.80 
8.07 
0.0 
0.0 



938.0 
550.0 
783.0 
0.0 
10.30 
12.60 
0.0 
0.0 
0.0 



M*/M 0.57 
K NM 211.7 

Usyra 43.5 



Nuclear Matter Characteristics 

0.60 0.67 0.54 

355.0 399.2 545.0 

36.1 43.9 35.0 



0.53 
626.3 



a The masses, the incompressibility Knm and the asymmetry energy a sym are in MeV, the coupling 
g3 in fm _1 . 
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TABLE II. RMF and constrained GCM results for the ground-state energies and mass rms 
radii and for the excitation energies AE n = (E n — Eq) of the first three monopole states calculated 
with the set NL1. 





Enerj 


;ies (MeV) 


Radii (fm) 


Excitation Enerj 


;ies (MeV) 


Nuclei 


RMF 


GCM 


RMF 


GCM 




AE 2 


AE 3 


16Q 


-127.24 


-127.46 


2.65 


2.65 


20.6 


38.9 


49.8 


40 Ca 


-342.48 


-342.58 


3.38 


3.38 


17.1 


29.9 


37.3 


90 Zr 


-784.90 


-784.99 


4.28 


4.28 


14.7 


29.1 


43.1 


208pk 


-1639.89 


-1640.07 


5.67 


5.67 


11.7 


23.3 


34.9 



TABLE III. Comparison of the GMR excitation energies (in MeV) obtained within the con- 



strained GCM calculations and the approximation (35) using the constrained incompressibil- 
ity (|34D with the nonrelativistic sum-rule approach obtained within nonrelativistic constrained 



Hartree-Fock (HF) calculations [32|. In the relativistic case sets NL1 and NL-SH are used, which 
have nearly the same nuclear matter incompressibility Knm as do the sets of Skyrme force param- 



eters SkM and SIII used in the nonrelativistic HF calculations [11| 





NL1: 


Knm = 


211.7MeV 


NL-SH: K NM 


= 354.95 MeV 




SkM: K NM = 


216.7 MeV 




SIII: K NM 


= 356.00 MeV 


Nuclei 


GCM 




Skyrme HF 


GCM 


Eq.® 


Skyrme HF 


16Q 


20.6 


20.9 


22.4 


25.3 


25.8 


26.6 


40 Ca 


17.1 


19.2 


20.2 


22.4 


23.9 


24.7 


90 Zr 


14.7 


16.3 


17.0 


20.16 


21.1 


21.2 


208p b 


11.7 


12.2 


12.9 


15.8 


16.1 


16.2 
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TABLE IV. Comparison of the RMF results from the parameter set NL-S1, with the HF+RPA 
calculations using density-dependent Skyrme interactions. Here we show the results only for the 
nuclei 208 Pb and 90 Zr, where the empirical data is reliable and rather well-established. 
Skyrme interactions+RPA GCM-RMF 





Dl 


Ska 


SIII 


NL-S1 


expt. 


Knm 


228 


263 


356 


296 




90 Zr 


18.5 


19.1 


22.1 


17.6 


17.0 ±0.5 


208p b 


14.4 


14.7 


17.2 


13.4 


13.5 ±0.3 
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FIGURES 

FIG. 1. The constrained energy Eq.(|33[) as a function of the rms radius R in Eq. (|3l|), for 90 Zr 
with the force NL-SH in the GCM calculations. The solid curve on the left is without taking the 
cut-off function into account. For the lower densities (higher R) the constrained energy has been 
extended (dot-dashed line) by inclusion of the cut-off function. The harmonic approximation to 
the curve is shown by the dotted parabola. The energies of the first three excited states (E\, E% 
and £3) along with the ground-state energy Eq are also shown. The energy E\ in the harmonic 
approximation (dotted line) differs only slightly from the actual value. 

FIG. 2. The ground-state vector and scalar densities for protons in the RMF (prmf) and in 
the GCM (poo) for 16 O. The effect of the ground-state correlations in the GCM densities are seen 
to be minimal. The GCM density for the first excited state of the GMR (pn) is also shown, where 
there is a considerable change in the densities in the interior. A slight change in the surface of the 
nucleus can also be seen. 

FIG. 3. The vector and the scalar transition density for the GMR in 208 Pb obtained in the 
relativistic GCM calculations with the force NL-SH. There is a conspicuous node in the densities 
of both the protons and neutrons. The change in the bulk of the vector density takes place at the 
expense of the change in the surface, thus conserving the total number of particles. 

FIG. 4. The constrained incompressibility Kq obtained in the RMF theory using various pa- 
rameter sets. Kc increases monotonically from NL1 to NL-SH for all nuclei. The values for HS 
show a slight dip, indicating a very large surface incompressibility. 

FIG. 5. The energy AEi of the GMR obtained with various relativistic Lagrangian sets in the 
GCM. The empirical values of the GMR in 208 Pb and 90 Zr have been shown at their average values 
by horizontal quadrangles. The widths of the quadrangles span the error bars in the empirical 
data. The empirical data encompass the corresponding GCM results at about K = 280 — 310 MeV 
and show a good agreement with the values obtained with the set NL-S1. 
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